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ABSTRACT 

I When a core collapse supernova occurs in a binary system, the surviving star as well 

(-H . as the compact remnant emerging from the supernova, may reach a substantial space 

Q^' velocity. With binary population synthesis modelling at solar and one fifth of solar 

metallicity, we predict the velocities of such runaway stars or binaries. We compile 
predictions for runaway OB stars, red supergiants and Wolf-Rayet stars, either isolated 
or with a compact companion. For those stars or binaries which undergo a second 
stellar explosion we compute their further evolution and the distance travelled until 
a Type II or Type Ibc supernova or a long or short gamma-ray burst occurs. We 
find our predicted population of OB runaway stars broadly matches the observed 
population of stars but, to match the fastest observed Wolf-Rayet runaway stars, we 
require that black holes receive an asymmetric kick upon formation. We find that 
at solar metallicity Type Ic supernova progenitors travel shorter distances than the 
progenitors of other supernova types because they are typically more massive and 
thus have shorter lifetimes. Those of Type IIP supernovae can fiy farthest about 48 pc 
on average at solar metallicity, with about 8 per cent of them reaching 100 pc. In 
' considering the consequences of assuming that the progenitors of long gamma-ray 

bursts are spun-up secondary stars that experience quasi-homogeneous evolution, we 
find that such evolution has a dramatic effect on the population of runaway Wolf- 
Rayet stars and that some 30 per cent of GRBs could occur a hundred parsecs or 
more from their initial positions. We also consider mergers of double compact object 
binaries consisting of neutron stars and/or black holes. We find the most common type 
of visible mergers are neutron star-black hole mergers that are roughly ten times more 
, common than neutron star-neutron star mergers. All compact mergers have a wide 

range of merger times from years to Gyrs and are predicted to occur three hundred 
times less often than supernovae in the Milky Way. We also find that there may be a 
population of low- velocity neutron stars that are ejected from a binary rather than by 
their own natal kick. These neutrons stars need to be included when the distribution 
of neutron star kicks is deduced from observations. 

Key words: gamma-rays: bursts - binaries: general - supernovae: general - stars: 
evolution - stars: Wolf-Rayet - stars: general 
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1 INTRODUCTION 

Runaway stars are isolated stars or binaries which 
have escaped from their parent clusters. A satisfactory 
set of observational characteristics that defines them 
is difficult to find. Often, it is assumed that a run- 
away star must have a space velocity of 30kms~^ 
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or greater. The most commonly o bserved massiv e run- 
away stars are OB runaways ([Blaauwl 196ll: Gied 
19871: IS tondll99l': ' Hoogerwerf. de Bruiine fc de Zeeuwll200ll : 



de Wit ct al. 2005) . Currently about 56 such stars are known 



in the Galaxy, (Hoogerwerf. de Bruiine fc de Zeeuwl 1200 ll ) 
with velocities up to 200 kms~^. There a re some cases know n 
with velocities greater than 500kms~^ (|Heber et al.ll2008l ). 
These high velocities and the stellar lifetimes of a few million 
years imply that such stars travel many parsecs from their 



2 John J. Eldridge 



initial positions before themselves exploding in core-collapse 
supernovae (SNe). 

There are two scenarios that create massive runaway 
stars. The first is dynamical ejection scenario (DES) where 
stars are ejected by close encounters in a dense cluster. 
The second is the binary supernova scenario (BSS) when 
one of the stars in a binary under goes a SN explosion 
and the system becomes unbound. The companion star 
then travels at roughly its pre-SN orbital velocity. It is 
thought that both scenarios contribute similar numbers of 
runaway stars. However evidence suggests that the BSS 
could be responsible for up to two thirds of obs erved run- 
aways (|Hoogerwerf. de Bruiine fc de Zeeuwll200ll '). Methods 
to discriminate between the two rely on the observational 
characteristics of the runaway star. Typically BSS runaways 
are expected to have a surface compositions that indicates 
they have experienced a binary interaction, while typical 
DES runaways might have main-sequence composition. A 
further complication is that some runaway stars, especially 
those with velocities in excess of about 200kms~^, may 
come from a binary that was ejected by the DES with one 
star being boos ted to high velocities after a BSS e jection 
as discussed by IPflamm- Altenburg fc Kroupal (I2OI0I ). Such 
binaries would be less common but could naturally explain 
the highest velocities observed for runaway stars. 

Interest in the final fat e of runaway stars ha s increased 
owing to the observations o f lHammer et all (120061). They ob- 
serve d that some long gamma-ray bursts ( Wooslev fc BloomI 
I2OO6I ) occur a few hundred parsecs away from the nearest 
star forming region. ICantiello et"ai] (|2007) suggest that this 
can be explained if the progenitor of the GRB was the sec- 
ondary in a binary system that was ejected via the BSS. 

Concurren tly the observations mad e by 

iFruchter et aT ('2006'), 'Kelly. Kirshner fc Pahrj (|2008r ) 
and Anderson & James (2008) also reveal details of how 
SNe are distributed in their host galaxies. They found that 
type Ic SNe and GRBs tend to be associated with the 
most luminous parts, and thus the sites of most intense 
star-formation, of galaxies while type II and type lb SNe 
tend to b e more evenly distributed t hroughout th ei r host 
galaxies. iLarsson et al.l (|2007h and iRaskin et al.1 (I2OO8I ) 
put this down to the difference in progenitor mass for the 
different SN types. The most massive stars explode earliest 
and thus closest to the sites of active star formation. While 
less massive stars have longer lifetimes and are not as closely 
associated with the most recent star formation. They did 
not consider the effect of runaway stars. It is important 
to understand how runaways change the distribution of 
SNe within a galaxy to remove any systematic error in 
determining the nature of different SN progenitors. 

In this paper we concern ourselves with describing the 
population of runaway stars from the BSS and their effect 
on the distribution of SNe within a galaxy. First we describe 
our method of simulating the runaways. Secondly we give 
our predictions for the space velocities of runaway stars, 
binaries and compact remnants and the distances our model 
runaways travel before exploding as SNe. Thirdly we outline 
how this affects how SNe are distributed with respect to 
their initial locations. Fourthly we perform a similar analysis 
for the merger of double compact object binaries owing to 
orbital decay by gravitational radiation. Finally we discuss 
our results and outline our conclusions. 



2 COMPUTATIONAL METHOD 

Our method is built upon the models and the popula- 
tion synthesis code described by lEldridge. Izzard fc Touti 
( 2008). We use the large number of detailed stellar evolution 
models that were calculate d with the Cambridge STARS 
evolution code, created by Egglet on (197 ^) and updated 
by vario us authors since (iPols et al.lil9951 : lEldridge fc Tout! 
20043; lEldridge. Izzard fc ToutI I2OO8I ). The population 
synthesis code then uses our detailed binary models to 
estimate various details of a binary population that are of 
interest and comparable with observations. For example, 
the relative numbers of different stellar types, the rela- 
tive rates of different SN types and, here, the expected 
velocities of runaway stars. Similar studies have used 
binary population synthesis to predict the runaway popu- 
lation before (Ivan Rensbergen. Vanbeveren fc de Loor3 
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.Vanbeveren. van Rensbergen fc Dc Looro 1998; iDrav et al.l 
l2005h . However this is the first time the effect of runaway 
stars on the distribution of SNe within galaxies has been 
considered. 

While a full description of our detailed models can 
be found in lEldridge. Izzard fc Tout! (|2008) we provide a 
brief overview. All the models employ our standard mass- 
loss prescription because it agrees best with vari ous ob- 
serva tions (Eldridge fcjbut 2004b; Eldridge, Izzard fc Tout! 
20081). For pre-WR mass loss, we use the rates of 



de Jager. Nieuwenhuiizen fc van der Huchtl (|l988D except 
for OB stars for which w e use the theoretical rates 
of IVink. de Koter fc LamersI (|200ll ). When the star be- 
comes a WR star Kurfa ce < 0-4, log f Teff/K ) > 

4.0], we use the rates of iNugis fc LamersI (|2000l ). We 
scale all rates with the standard factor (Z/Z ^)"'^ 
(|Kudritzki. Pauldrach fc Puis! 1 19871: iHeger et al.ll2003l ). ex 



cept for the rates of Vink. de Koter fc Lamer"^ 200 ih which 
include their own metallicity scaling. 

We have modified our stellar evolution code to model 
binary evolution. The details of our binary interaction algo- 
rit hm arc relatively simple com pared to the scheme outlined 
bv iHurlev. Tout fc PoS (|2002l ). We used their scheme as a 
basis but we changed some details which cannot be directly 
applied to our detailed stellar evolution calculation. We also 
make a number of assumptions to keep our code relatively 
simple. Our aim was to investigate the effect of enhanced 
mass loss due to binary interactions on stellar lifetimes and 
populations. Therefore, we concentrated on this rather than 
every possible physical process which would add more un- 
certainty to our model. We also make assumptions in calcu- 
lating our synthetic population to avoid calculating a large 
number of models. For example, we do not model the accre- 
tion on to the secondary in the detailed code. We take the 
final mass of the secondary at the end of the primary evolu- 
tion as the initial mass of the secondary when we create our 
detailed secondary model. This avoids calculating 10 times 
more secondary models than primary models. 

We always define the primary as the initially more mas- 
sive star and we only evolve one star at a time with our 
detailed code. When we evolve the primary in detail, it has 
a shorter evolutionary time-scale than the secondary which 
remains on the main sequence until after the primary com- 
pletes its evolution and so we can determine the state of 
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the secondary with the sin gle stellar evolution equations of 
iHurlev. Pols fc Tout! (|2000l ). When we evolve the secondary 
in detail, we assume that its companion is the compact rem- 
nant of the primary (a white dwarf, neutron star or black 
hole) and treat this as a point mass. 

If Roche lobe overflow occurs mass lost from the pri- 
mary is transferred to the secondary but not all is neces- 
sarily accreted. Accretion causes the star to expand owing 
to increased total mass and therefore an increased energy 
production rate if M2 A/2/TKH, where tkh is the thermal, 
or Kelvin-Helmholtz, time-scale. We assume that the star's 
maximum accretion rate is determined by its current mass 
and its thermal timescale. We define a max;imum accretion 
rate for a star such that Af2,max = M2/tkh- If the accretion 
rate is greater than this, then any additional mass and its 
orbital angular momentum are lost from the system. In gen- 
eral, stars with lower masses have longer thermal time-scales 
than more massive stars. Efficient transfer is only possible 
if the two stars are of nearly equal mass so the thermal 
time-scales are similar. This is an approximate treatment 
but provides a similar result to th e more complex model of 
IPetrovic. Langer fc van der Huchtl ((2005 ) who included ro- 
tation and found that it led to inefficient mass transfer. For 
compact companions, we d erive the maximum accre tion rate 
from the Eddington limit l| Cameron fc Mocklll967l ). 

In summary our stellar models are a set of around 
15,000 detailed stellar evolution models of single stars, 
primary stars, secondary s tars and merged systems 
l|Eldridge. Izzard fc Tou^l2008l '). We also include a new se- 
ries of models to account for quasi-homogeneous evolution. 
These are discussed below. For this work we made several 
refinements to our population synthesis code to predict run- 
away velocities and account for the possibility of SN reversal. 
That is when the initially lower-mass secondary star under- 
goes its SN first after mass accretion. 

2.1 Population synthesis 

Our population synthesis calculations are b uilt upon those 
described bv jEldridge. Izzard fc Tou t! (|2008l ) with some im- 
provements. Here we focus on the evolution after the first 
SN in the binary. First we estimate the lifetimes of the pri- 
mary and the secondary star. Our secondary models only 
have a limited companion mass range but this provides a 
reasonable estimate because the lifetime varies little with 
different amounts of mass loss. We compare the lifetimes of 
the two stars to check whether the SN order is reversed. We 
find that it is more likely at higher metallicity but only oc- 
casional, occurring in 7 per cent of our binary systems. If 
the secondary lifetime is shorter than the primary lifetime 
we perform a similar calculation as described below but now 
the secondary explodes first and may eject the primary star. 
This gives rise to more stripped SNe occurring away from 
their initial position. 

If a star has a carbon/oxygen core mass greater than 
1.38M© and the final mass of the star is greater than 2Mq 
we assume it explodes in a SN. We select the SN type as de- 
scribed below. We estimate what type of compact remnant 
a stellar model will produc e by using the method outline 
in lEldridee fc Tout! (|2004bh . We assume that first a neu- 
tron star is formed at the centre of the star after core col- 
lapse of mass Mch ~ IAMq. This produces about 10'*® J 



of energy from the release of gravitational binding energy 
in neutron star formation. We then assume a hundredth 
of this energy is transferred into the envelope by some un- 
known mechanism. The current suggestion is the transfer 
occurs via neutrinos released from forming the protoneu- 
tron star that are thermalized within the envelope or dense 
outer parts of the core. We integrate the binding energy 
of the star from the surface towards the centre until we 
reach 10** J. The envelope outside this region is ejected 
with the remaining amount forming the remnant. If we have 
Mrem > SM© , it is & black hole and we set Mbh = Mrcm- 
Otherwise we have a neutron star with mass IAMq. We de- 
termine the fate of the binary if a neutron st ar is formed 
by th e work of iTauris fc Takend (figQa ) and iTauris et al.l 
()l999l ) with the latest determin ation for the kick v elocity 
distribution from observations of iHobbs et al.l l|2005h . If the 
system is unbound then t he velocities of both sta rs are cal- 
culated by the method of ITauris fc TakensI (|l998l ) that con- 
siders every relevant factor. However we neglect the super- 
nova impact on the companion star. If the system remains 
bound then the veloc ity of both components is the resul- 
tant system velocity (|Brandt. Pod siadlowski & SigurdssonI 
I1995I ). If the remnant is a black hole, we assume that it 
receives a similar kick. Because the masses of black holes 
are greater t han those of neutro n stars we use the kick dis- 
tribution of iHobbs et al.l (|2005l ') but as a momentum dis- 
tribution. We calculated the black hole kick velocity, vbu 
from vbh = fNs(l -4^/0 /Mbh). Where uns is a kick veloc- 
ity selected at random from the neutron star kick velocity 
distribution and Mbh is the mass of the black hole. Black 
holes are not normally considered to have significant kicks 
although their importa nce has been invest i gated by others 
dBrandt. Po"dsra dlowski fc Sigurdssonlll995l: iKaloeeral Il999l : 
iMirabel et a'Lll200a . l20o3 : IVoss fc Taurisll2003l ). Wefind that 
we must include them or we obtain no runaway WR stars 
with initial masses greater than 30Mq. 

The resulting velocities of runaway stars and runaway 
binaries are recorded. We also record the SN type and its 
location at the initial position of the binary. We weight 
each event by a Saltpeter IMF and assume the distribu- 
tions of initial mass ratios and the logarithm of separa- 
tions are fiat. Our models have initial separations that 
take values between 1 < logiQ(a/7?Q) < 4 in steps of 
0.25 dex. The mass ratio takes values of q — 0.3, 0.5, 0.7 
and 0.9. We do not use the q — OA models employed by 



lEldridge. Izzard fc ToutI (|2008h because there is growing ev- 
idence that t he mass ratio is skewed to larger values in mas- 
sive binaries |Garmanv et al.lll980l: iPinsonneault fc Stanekl 



I2OO6I : iKobulnickv fc Frveill2007l : iKiminki et al.ll2009l ). 

We next consider the fate of the secondary star. If it 
accreted material from the primary we assume the star is 
rejuvenated and use its post mass transfer mass as its new 
initial mass. If a secondary accreted more than 5 per cent 
of its initial mass when it was a main-sequence star we as- 
sume it has been spun up and rotationally induced mixing 
mixes fresh hydrogen into the core and it is rejuvenated 
and take it to be a zero-age main-sequence star. We do 
this at all metallicities. This extends the lifetime of some 
of our runaway stars. This is a similar assumption as used 
by IVanbeveren. vaii Rensbergen fc De Loor j (|l998D . If the 
binary was unbound by the first SN we use a single star 
model to determine the result of its evolution. If the system 
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remains bound then we use our secondary models. If the 
secondary experiences a SN it does so after it has travelled 
away from the location of the primary SN. We determine the 
distance travelled by considering in detail the geometry of 
situation using the runaway velocity relative to the original 
plane, the orientation of the binary to the line of site, the 
phase of the stars in the binary and the time that has passed 
since the first SN. We then record the location of the SN. 
Finally we determine the final outcome of the binary and 
whether the system is unbound or a double compact object. 
We record the velocities of the single and binary compact 
object. 

We perform the above analysis over the full range of our 
primary models. Because we select our neutron star kicks at 
random, we repeat our calculations a large number of times 
to cover the full range of possible outcomes after the first 
SN. 

2.2 Quasi- homogeneous evolution 

We include one new evolutionary path, not in the stan- 
dard picture of binary evolution for our secondary stars. 
If a secondary accreted more than 5 per cent of its initial 
mass as discussed above we assume it has been spun up 
and rotationally induced mixing mixes fresh hydrogen into 
the core and it is rejuvenated and take it to be a zero-age 
main-sequence star. If the secondary's initial metal mass 
fraction is less than or equal to 0.004 and it has a mass 
after accretion of more than WMq, we assume it contin- 
ues to evolve fully mixed during its entire main-sequence 
lifetime. This is referred to as quasi-homogeneous evolution 
(QHE). It is the result of rapid rotation due to the accre- 
tion of materia l from the primary s tar an d is described by 
[Maedcr ('1987'), 'Mevnet & Maeder (200?), 'Yoon & Langei) 
( 2005), Yoon, Langer fc Norman (^2006) and .Cantiello ot al. 
(|2007l ). The star does not spin down as at lower metallicity 
stellar winds are weaker so less angular momentum can be 
lost. 

To include we use simple models in which we assume 
the stars are fully mixed during their hydrogen burning evo- 
lution. Once hydrogen burning ends this extra mixing is 
turned off. We use the models whether the binary was bound 
or unbound after the primary SN. We find, a.t Z — 0.004, 
QHE increases the percentage of SNe that are type Ib/c 
from 20 to 26 per cent. We find that three per cent of all 
stars that experience QHE give rise to a long-GRB as we 
also require a final CO core mass of greater than 7Mq for 
a long-GRB to occur. More metal rich stars do not experi- 
ence QHE as we assume the stars rapidly spin-down as the 
stronger stellar winds take away angular momentum more 
rapidly. 

Our requirement for QHE to occur may seem quite re- 
laxed. [Canticllo ct al. ( 2007,) restricted their study to a Case 
B, post-main-sequence mass-transfer system rather than 
Case A system. This was to avoid orbital synchronization 
slowing the rotation of the secondary. Therefore we may be 
overestimating the number of long-GRBs via QHE. However 
most of the QHE systems in our models are Case B and also 
systems with high mass ratios which have similar thermal 
timescales so mass-transfer is most efficient. Therefore the 
opportunity of tidal synchronization is low. Further more 
such synchronization is unlikely as the stars have radiative 



Table 1. Tlie relative population of mergers from our binary 
population at different metallicities and with different minimum 
masses for the binary systems. 



Mass Range 




During 


Post-MS 


No Merger 




Z 


MS 






5 -> 120 


0.004 


0.029 


0.001 


0.970 




0.020 


0.036 


0.007 


0.957 


10 120 


0.004 


0.070 


0.002 


0.928 




0.020 


0.088 


0.016 


0.896 


20 -> 120 


0.004 


0.106 


0.005 


0.889 




0.020 


0.129 


0.004 


0.867 



envelopes during much of their evolution. Thus any tidal 
forces are likely to be weak. 

In addition our 5 per cent increase in mass required for 
QHE may also seem arbitrary. Only about 10 percent of the 
mass of a main sequence star needs to be accreted to bring 
its spin from zero to critical. However QHE does not need 
critical rotation and the initial rotation rate of these stars 
is not zero. Therefore 5 per cent is a reasonable limit. A full 
treatment would require a more detailed model of rotation 
within the star and the binary. However the inclusion of this 
would severely limit the stability of the code and thus make 
it currently difficult to create the large number of models 
required for population synthesis. 

2.3 Non-degenerate mergers 

In Table [T] we list the relative fraction of binaries that merge 
before the first SN occurs. We find that 7 to 10 per cent of 
binaries with primaries massive enough to explode in a SN 
experience a merger. Therefore the majority of binaries pro- 
vide two possible supernova progenitors and contribute to 
the runaway stars. However even stars not massive enough 
to provide a SN can merge and in some cases the result- 
ing star is massive enough to produce a SN. We find more 
mergers in the higher mass binaries. 

2.4 Determining SN types 

Determining what SN type a star produces is a difficult 
task. Here we use the relative SN rates of ISmartt el al.l 
(|2009l) to estimate the parameters required for the different 
types. The rates are based on a volume limited sample of 
92 core-collapse SNe within 28Mpc over 10.5 yrs. The range 
of metallicities for the SNe tend to be an even mix of LMC 
and solar like {Z — 0.008 and 0.020 respectively). To match 
the relative rates, we first select out stars that we expect 
to explode in supernovae. We first require that core carbon 
burning has occurred and that the total mass of the star 
is greater than 2Mq and that the carbon-oxygen core mass 
is greater than l.SSM©. Next we vary the total masses of 
hydrogen to helium and their relative amounts required for 
the different SN ty pes. To reproduce the observed rates of 
ISmartt el all (|2009l ) we find that we must include one single 
star for every binary system included in the code. Otherwise 
we cannot reproduce enough type IIP SNe while keeping the 
maximum type IIP progenitor mass below 2OM0 . This is the 
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99 per cent confidence limit o n the maximum ma ss of type 
IIP progenitors determined bv lSmartt efall (|2009l ') from ob- 
servations. This in effect increases the number of very wide 
non-interacting binaries in our sample. We summarize the 
parameters required for each SN type in Table (2] If a type 
II SN does not produce a IIP event we include it in a broad 
cl ass of non-IIP SNe. These events are rare and the statistics 
of lSmartt el"aLl (|2009l ) are not good enough to determine the 
parameters required for the sub-types of lib, Iln and IIL. 

For the hydrogen-poor SN-types we differentiate be- 
tween type lb and Ic, helium rich and poor SN respec- 
tively, by considering the amount of helium in the ejecta. 
Rather than u sing the total mass of helium in the ejecta as 
other authors (WcUstci n fc Langeilll999l:|Pols fc Dewill2002l : 
iGeorev et alll2009. : .Yoon. Wooslev fc Langerll2010l ) we con- 
sider the fraction of helium in the ejecta to differentiate be- 
tween lb and Ic SNe. We prefer this method because oth- 
erwise some progenitors with pure helium ejecta of only a 
fraction of a Solar mass would be identified as Ic rather 
than lb. This means our results differ from previous authors. 
The main difference is that while other authors find two dis- 
tinct classes of low and high mass type Ic progenitors with 
type lb progenitors of mainly low masses, we find more of a 
overlapping c ontin uum of type Ic and lb progenitor masses. 
iDrout et al.l (|201ll ) suggested this might be the case from 
a systematic study of lightcurves of type Ib/c SNe. At low 
metallicity we note more progenitors explode as type lb SNe. 
The latter is due to weaker WR winds removing less helium 
before core-collapse. 

The observed relative rates of supernovae from 
ISmartt el al.l l|2009l ) are shown along with our synthetic pop- 
ulation rates from an equal mix of solar and LMC metallicity 
in Table [3] The predicted ratios and mean initial masses for 
the different supernova types are given in Section 3.3. We 
find that the mean initial mass of the progenitor stars in- 
creases through the SN types from IIP being the least mas- 
sive to Ic progenitors the most massive. The mass ranges 
of the different SN types are shown in Fig. [1] Type IIP SN 
progenitors are not only here the lowest mass but also the 
narrowest range of masses. The remaining SN types have 
larger ranges of both initial and final mass. 

As mentioned above we also record SNe that 
may have an associated long-Gamma-Ray Burst (GRB) 
l| Wooslev fc Bloomll2006l ). While we do not consider rota- 
tion in our models, to produce a GRB we require it must 
have experienced QHE as described above, that the CO core 
mass must be greater than IMq and that the star's rota- 
tion axis is within 10 degrees of the line of sight. The direc- 
tion of the rotation axis is taken to be random. We assume 
that the rotation axis is perpendicular to the orbital plane. 
This is important to consider when we calculate the appar- 
ent distance travelled by progenitors before they explode as 
discussed below. Due to our assumptions about QHE, long- 
GRBs only occur at a metallicity of Z = 0.004. One in ten 
QHE stars leads to a long-GRB and there is one long GRB 
for every 175 type Ib/c SNe or 685 of all core-collapse SNe. 
Also we find that the mean initial mass is lower than the 
mean mass for type lb and Ic SNe. This is because the stars 
must accrete mass to be spun up and this biases the progen- 
itors towards lower masses. 



Table 2. The required parameters for a star to 
specific SN type. 



dve rise to a 





Final 


CO core 








SN 


Mass 


Mass 


M(H) 


M(H) 


M(He) 


Type 


/Mq 


/Mq 


/Af(Hc) 


/Mq 


/M(ejecta) 


IIP 


> 2 


> 1.38 


> 1.1 


> 0.05 




II 


> 2 


> 1.38 


< 1.1 


> 0.05 




lb 


> 2 


> 1.38 




sC 0.05 


^ 0.58 


Ic 


> 2 


> 1.38 




s; 0.05 


< 0.58 



Table 3. Relative fraction s of different SN types from observa- 
tions of lSmartt el al.l I^Oo3) and output from our population syn- 
thesis with a mix of single stars and binaries at a mix of metal- 
licities. 



z 


IIP 


non-IIP 


lb 


Ic 


0.008 & 0.020 


0.588 


0.122 


0.095 


0.195 


Smartt et al. 


0.587 


0.120 


0.098 


0.195 



3 RESULTS 

While our code can run at any metallicity we restrict our- 
selves to just two, solar metallicity with a metal mass frac- 
tion Z = 0.02 and SMC-like metallicity with Z = 0.004. 
This enables us to decouple the effects of metallicity on our 
results. 



3.1 Single-star Runaways 

The main output of our code is the distribution of runaway 
velocities. Observations of runaways have been used to infer 
that around 10-30 per cent of O stars and 2-10 per cent 



Table 4. The fractions of OB stars, red supergiants and Wolf- 
Rayet stars that are runaways in our synthetic population. We 
list the fractions from two selection criteria for runaways. First 
any star that has a peculiar velocity above 5km s~^ and a sec- 
ond with a higher velocity requirement of 30kms~^. The former, 
lower limit, would provide results reco rds runaways in a method 
similar to the observational method of IStonj l ll99ll l. We list two 
populations, one from our standard population including single 
stars and binaries, and a second including some of the single stars 
as DES runaways, we assume half the single stars are runaways 
for half their O star lifetime. 
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Figure 1. The probability distribution of final (pre-explosion) and initial masses for different type SN progenitors from our population 
including single stars and binaries. The distribution of initial masses are shown by the dashed lines and the final masses are shown by 
the solid lines. The models on the left have a metallicity mass fraction o{ Z = 0.004 and the models on the right have a metallicity mass 
fraction of Z = 0.020. 



z = 0.020 




Velocity/km s ' 

Figure 2. The cumulative frequency of stellar runaway 
velocities. The thick black line is th e observations from 
iHoogerwerf. de Bruiine &: de Zeeuw* f200l"), the thick dashed line 
is for the same data but with the known binary stars removed, 
the thin dashed line is the line predicted by our models for all 
runaway stars, the blue line for all OB stars, the red line for 
all red supergiants and the green line for Wolf-Rayet stars. The 
predictions are for solar metallicity. 

of B stars are observed to be runaways (lGieslll987l : IStond 
Il99ll : Ide Wit et all l2005l : IZinnecker fc Yorkd l2007i v These 
observed fractions are very uncertain as the number of well 
observed runaways is relatively small and the definition of 
what a runaway is varies between different studies. For ex- 



ample the above observed runaway fractions vary b ecaus e 
the runaway defi nition varies between authors lciesl l|l987l ). 
IStond (liggih and Ide Wit et al.l lj2005h state that a runaway 
must have a velocity greater than 30km s~^ or be a large 
distance above the Galactic plane. Also in some cases is it 
not clear if the percentage given is the number runaways 
in the field O star population or the number of runaways 
relative to the entire O star population. 

The recent study of Ide Wit et al.1 l|2005h suggest the 50 
per cent of field O stars are runaways. They also state that 
since 70 per cent of O stars are in clusters. These facts sug- 
gest that possibly 15 per cent of O stars are runaways. They 
identify runaways by peculiar velocities above 40km s"'^, a 
high distance from the galactic plane or proximity to young 
clusters. These mixed definitions make it difficult to identify 
a velocity cut-off to use with our model population to pre- 
dict the number of runaways. For example even a runaway 
with a velocity of 5km s~^ would travel approximately 10 
pc in 1 Myr and w ould appear as a runaway according to 
Ide Wit et all (|2005l ). Therefore we estimate a range of the 
fraction of runaways from our models. We use the high ve- 
locity of 30km s~^ to give a strict minimum runaway fraction 
estimate and a lower velocity limit of 5km s~^ to provide a 
more relaxed maximum runaway fraction estimate. 

In Table |4j we list the fraction of the predicted popula- 
tions that are runaways. The numbers are calculated weight- 
ing each model by its lifetime in the phase of evolution and 
with a Salpeter IMF. For comparison we list the results for 
O stars, B stars, red supergiants (RSGs) and WR stars. B 
stars are taken to be stars with \og^f^{Tes / K) ^ 4.18 (B5 and 
earlier only includes stars with initial masses greater than 
5M0, the minimum mass of our binary models) and a sur- 



Runaway stars as progenitors of SNe and GRBs 7 



OB stars z = 0.020 



WR stars z = 0.020 



100 - 




30 40 60 80 100 150 200 250 

Velocity/km s"' 

OB stars z = 0.004 



100 - 



1 1 


1 1 1 






























\ N \ ^ \ 












^^^^ 1 ^\ 1 ^\ 











30 40 60 80 100 150 200 250 

Velocity/km s"^ 



100 



10 





1 1 1 










- ^^^^^^^^^^^ 










y/ 


1 III II 



30 40 



60 80 100 150 200 250 
Velocity/km s"^ 

WR stars z = 0.004 



100 - 




30 40 60 80 100 150 200 250 

Veloclty/km s"' 



Figure 3. Contour plots of the probability for runaway star to have a certain mass and velocity. The units of the contours are probability 
per 10kms~^ and per 0.1 in \o^-^q(M IMq). The mass referred to is the actual mass so includes any mass gained b y a runaway during 
mass transfer. The red points are the observations of OB runaways from lHoogerwerf. de Bruiine &: de Zeeuwl ll200lh with stellar masses 
derived from fitting the spectral energy distribution. The top two panels are for solar metallicity and the bottom two panels are for 
SMC-like metallicity. The left hand panels are for OB stars while the right hand panels are for WR stars. The mass axis is either the 
initial mass of the runaway star or its mass after accretion due to any binary interaction. A similar figure showing the number of stars 
below the velocity limit for a star to be observed as a runaway, 30km s~^, is included in the appendix (Fig. [T2J. 



face hydrogen mass fraction greater than 0.4, O stars have 
hotter effective temperatures of logjQ(Tcff /Jf) > 4.48 and 
4.52 at Z = 0.020 and 0.004 respectively. RSGs are defined 
by logi(,(rcff/A') ^ 3.66 and Iogio(L/LQ) > 4.9 and WR 
stars are taken to have a surface hydrogen mass fraction 
less than 0.4, log^Q(Tcff/Ar) > 4 and the same luminosity 
limit as the RSGs. 

In calculating these numbers we consider count each O, 
B, RSG or WR star including both the primary or secondary 
in the binary. It is unclear if in the observed samples an O 
star binary is counted as one or two O stars. We find 0.5 to 
2.2 per cent of O stars are runaways and 2.2 to 7.1 per cent 
of B stars are runaways at Solar metallicity. 

We note that uncertainties in our model may mean the 
O runaway fraction is an underestimate. For example, the 
absence of rotational mixing in our stellar evolution mod- 
els. Rotational mixing mixes fresh hydrogen into the core of 



a main-sequence star and this would increase our number 
of O star runaways by increasing their lifetimes. The evi- 
dence that rotation has an importance in the number of O 
star runaways if we include QHE, the most extreme form 
of rotational mixing, at solar metallicity the number of O 
star runaways increases to 2.2 to 16 per cent. We note the 
dramatic increase in number of O runaways is due to stars, 
of around 15Mq that have 20 Myrs O star lifetimes with 
QHE but are not O stars without it. Rotation will have a 
similar but less dramatic effect at Solar metallicity. The ef- 
fect of rotation on a stars lifetime is only to increase it by 
approximately 10 per cent. Including a detailed model of 
accretion on to the secondary, especially following rotation 
would allow more accurate inclusion of the effects of reju- 
venation on the secondary stars that accrete material from 
their primaries. 

The number of B star runaways is similar to that sug- 
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gested from observations. We seem to under-predict the 
number of O star runaways if we compare to the typical num- 
be rs quoted for the number of runaway O stars. The work 
of lde Wit etin |[2005il used various O star catalogues in- 
cluding that of Maz-Apellniz et al.l (|2004| ) . We use the same 
catalogue to estimate from the catalogue the fraction of O 
star runaways. It is important to note that we count O star 
binaries in the catalogue as 2 O stars. We find that 76 ± 5 
per cent of O stars are found in clusters with 6±1 per cent of 
O stars being runaways. The list of runaways includes some 
binaries, we assume these are DES in origin and find that 
there are 4 ± 1 per cent of single O star runaways. Consid- 
ering the uncertainties in classification of runaways and the 
approximations within our population model our predicted 
runaway population is comparable to the runawa y popula- 
tion in the catalogue o f iMaz-Apellniz et al.l |200i). 

Further complicating matters is that we have not con- 
sidered in our model the contribution of DES runaways. We 
have made a toy-population assuming that our single star 
population provide a similar number of DES runaways to the 
total number of BSS, approximately 30 percent of the binary 
produce BSS. We then assume that they are runaways for 
half their O star lifetime and have a similar velocity distri- 
bution to the BSS runaways. We see in TableUthese modest 
assumptions increase the number of O runaways to between 
1.5 to 6.6 and the number of B runaways to between 2.8 
and 9 percent. These numbers are also comparable to the 
observed numbers. 

After considering the population of runaways we now 
consider the velocity distribution and parameters of indi- 
vidual runaways. The b est observed sample of OB runaways 
to da te is that from iHoogerwerf. de Bruiine &c de Zeeuwl 
who list the velocities of 56 runaway stars. 
IHoogerwerf. de Bruiine fc de Zeeuw* ('2001') give mass es for 
some of these estimated by Schm idt-Kalo r, (|l982l ') and 
IVanbeveren. De Loore fc Van RensbergenI (|l998l V By us- 
ing available UBVJHK p hotometry we have empl oyed 
the methods outlined in lEldridge fc Relanq (|2010l ) to 
estimate stellar masses with Cambridge STARS mod- 
els. Our masses broadly agree with those used by 
IHoogerwerf. de Bruiine fc de Zeeuwl (|200ll ) and are used in 
Fig. (31 The mass range we find is between 5 to 70Mq with 
half the stars having masses less than IOMq . 

In Fig. [2] we plot the cumulative frequency of our run- 
away population velocities. In general there is agreement be- 
tween the observed runaway velocity distribution and that 
predicted by our models. We find that the assumed initial 
binary separation distribution largely determines the shape 
of the runaway velocities. Here we use a distribution that is 
flat in logj^Q a. If we limit ourselves to only close binaries we 
over predict the number of fast runaways, while considering 
the widest binaries leads to only the slowest runaways. It ap- 
pears in the flgure that we are under-predicting the number 
of fast runaways. 

Using our estimated initial masses from the runaways 
we can extend Fig. [2] along another axis as in Figure [3] of 
effective initial mass of the runaway star. This is to take 
account of the increase in its initial mass due to accretion of 
mass from its primary compani on. The general shape o f our 
predictions is similar to that of lPortegies ZwartI (|2000l ) who 
presents a similar flgure. In general it is clear that less mas- 
sive stars are more likely to achieve higher runaway veloci- 



ties. The observed runaway stars agree with this prediction 
however there are several outliers with high masses and mod- 
erate velocities. Observations of these stars indicate they 
also have large rotational velocities , v sin i ~ 200kms~^ 
(|Hoogerwerf. de Bruiine fc de Zeeuwl 120011 ). and therefore 
their masses from SED fltting should be considered upper 
limits because they were derived with non-rotating single 
star models. Rotational mixing can increase the luminosity 
of a star and the main-sequence lifetime. This means their 
mass would be overestimated and make them more likely 
to be observed. Alternatively, they could be DES runaways, 
runaways that experience both DES and BSS or our binary 
models are producing mergers in very close binaries rather 
than producing fast BSS runaways from close binaries. 

We must be wary of the selection effects of the ob- 
served sample. It is easier to observed higher velocity and 
high mass runaway stars. These selection effects should be 
more carefully understood before any detailed comparison 
between observations and theory. For example most run- 
aways are likely to have masses less than IOMq with ve- 
locities below lOOkms"^. However such stars and velocities 
are difficult to detect within the Hippa r cos ca talogue used 
by IHoogerwerf. de Bruiine fc de Zeeuwl (|200ll ). Meanwhile 
more massive, more luminous stars are easier to detect, es- 
pecially if they have large proper motions. Quantifying these 
selection effects is difficult. We may need to wait for the Gaia 
satellite before obtaining a less biased catalogue of runaway 
stars. Because the observed population is uncertain we do 
not attempt further tuning of our initial binary distribution. 

We note that at both metallicities the fraction of sec- 
ondary stars that have accreted some material during mass 
transfer is 50 per cent. In some cases the amount accreted is 
very little, only 17 per cent of secondary stars accrete more 
than 5 per cent their original mass due to assumptions on the 
maximum accretion rate for the secondary. Because of this 
we estimate that no more than 50 per cent of BSS runaways 
should be observed to have high rotation velocities and/or 
enriched compositions. Therefore, absence of these observ- 
able features should not be inferred to mean the runaways 
did not come from BSS. 

Finally we again note that our predictions are for the 
BSS only. However, the main detail that determines the ve- 
locity distribution in both cases is the initial binary separa- 
tion distribution assumed. We therefore suggest that the ex- 
pected velocity distribution from the DES should be similar 
to that of the BSS. Stars that are ejected by DES normally 
follow the interaction of a binary with anoth er binary or a 
single star. iGvaramadze fc Gualand"ri3 (|2010l ) have recently 
investigated the runaway velocities possible from three body 
interactions. They found that the mean velocities are similar 
to those we show here in Figure [2| However they also find 
the interactions can lead to velocities in excess of 100 kms~^ 
for BOM© stars in 10 per cent of triple body encounters. This 
adds weight to the suggestion that runaways more massive 
than 4OM0 with runaway velocities faster than 40kms~^ 
are in fact DES or a DES-I-BSS combination. 



3.1.1 Wolf-Rayet runaways 

Fig. O shows our results for WR stars. The sharp cut off at 
low masses in these panels is due to our luminosity limit 
that WR stars must have logjo(I//Z/0) ^ 4.9 to be included. 
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The difference between the two metalhcity plots at higher 
masses is due to the differing lif e times of WR stars at the two 
metallicities. lEldridee fc VinkI (|2006') found that, at SMC- 
like metallicity, the lifetimes of WR stars vary between 3 
to 4x10^ yrs while at solar metallicity the lifetimes are 4 
to 8x10^ yrs. Furthermore QlfE also skews this plot be- 
cause such stars become WR stars while still on the main 
sequence and so the plot begins to resemble the OB star 
panel as the time-scales are millions of years rather than a 
few hundred thousand years. If we do not include QHE then 
the distribution becomes the same as for the solar metallicity 
plot. Thus a large population of runaway WR stars in low- 
metallicity galaxies with an average velocity above 50 km s~^ 
would be circumstantial evidence for the occurrence of 
QflE in nature. Surveys for runaways of the SMC have 
be u ndertake n but have not yet provided detailed statistics 
iGva ramadz e. Pflamm-Altenburg fc Kroupall2010l ') . 

The fraction of WR and RSG stars that are runaways 
is larger than the fractions of OB runaways (Table |4|. This 
is because BSS runaway stars spend a large fraction of their 
main-sequence life stationary, while they spend most or all 
of the lives as RSGs or WR stars as runaways. This makes 
their relative number appear greater. At SMC metallicity, 
the WR runaway fraction is particularly high due to QHE 
scenario, which produces more WR stars from secondary 
stars. 

The velocities of Wolf-Rayet stars differ strongly be- 
tween the two metallicities we study. At solar metallicity we 
find WR runaways have velocities less than about 150 km s~^ 
but runaways with velocities above 80kms~^ are rare. If we 
do not included black hole kicks we only find WR runaways 
with initial masses less than SOM© and 0.4 per cent of WR 
stars have velocities above 30kms^^. This is because sec- 
ondaries that become WR stars owing to stellar wind mass 
loss have primaries that form black holes at core-collapse 
and remain bound because very little mass is ejected in such 
supernovae. We find that WR stars eject relatively little of 
their final mass, between 2 to SM©. With such low ejecta 
masses without a strong kick it is difficult to unbind the 
binary. It is possible that, due to the sparseness of our bi- 
nary grid, we may be missing WR runaways from binaries 
in which both components have initial masses in the range 
20 to 25M0. We might also assume too low an initial mass 
limit for a black hole to be formed in core collapse. Binaries 
with component masses between 20 and 25M(7) may give ris e 
to high velocity runaways as shown by iDrav et ahl (|2005h , 
fig. 2(c), which refers to a similar population synthesis with 
no black hole kicks. In these systems the binary mass ratio 
at core collapse is more even so the secondary may achieve 
a higher velocity. Such systems are still rare and dominated 
by binary WR stars remaining in their binary systems. 

Our re sults br oadly agree with similar theoretical pre- 
dictions of lDrav"e t al. (200^). It is difficult to perform a 
quantitative comparison but in their fig. 2(b), which de- 
scribes their model most similar to ours as it includes black 
hole kicks, we find broadly similar maximum velocities for 
WR stars of up to around 120 km s ~^. However most WR 
stars have velocities below 80 km s" - MOrav et all (|2005l ) also 
find a restricted number of WR runaways when they assume 
that black holes do not have kicks. 

We note that iDrav et all (|2005l '). in fig. 2(a) predict a 
large number of runaway WR stars. In this model they ig- 



nored stellar-wind mass loss. So more of the secondary stars 
accrete enough material to become WR stars as single stars, 
when the stellar wind mass-loss rates are reapplied after the 
first SN. Such an arrangement is unlikely to occur in na- 
ture because, if the metallicity were low enough to reduce 
the main-sequence mass-loss rates, the minimum mass for a 
WR star would also increase. The only way for WR stars to 
occur in such a situation would be via QHE. 

Another reason for mostly slow WR runaways in our 
synthetic population is demonstrated by the example that 
5OM0 primary stars at solar metallicity form neutron stars 
in core collapse. These binaries, even with large neutron star 
kicks, only obtain slow runaway velocities because the sec- 
ondaries are much more massive than their companions at 
core-collapse and so their orbital velocities are smaller. 

Finally we may also predict more WR runaways if we 
were to relax our luminosity constraint for a WR star. This 
number is highly uncertain and if it is set too low we would 
overpredict the number of WN stars. This would suggest 
that the higher velocity a runaway WR star has, the lower 
its initial mass. 

Observa t ions of runaway O and WR stars by 
I Moffat et al.l (| 19981 ) imply that there are a similar fraction 
of runaways for both types of star. They also suggest that 
mean kinetic ages based on displacement and motion away 
from the Galactic plane tend to slightly favour the DES 
over BSS. There are a few WR stars with higher veloci- 
ties such as WR 124 with a space velocity of 180kms~^ 
(Ivan der Sluvs fc Lamerd [2003 1 . Given our results (Fig. [3]), 
it seems unlikely that WR 124 was ejected through the 
BSS. Conceivably, it may be the res ult of both DES an d 
BSS combined. Or, as suggested by iMoffat et al] (|l998l) , 
the more massive a runaway the more its probability of 
coming from DES incre ases. Recent radio observation by 
iDzib fc Rodrigue3 (|200!jh give velocities of about 30 kms~^ 
for a few WR stars, in agreement with our predictions. 
However they also give a runaway velocity of more than 
300 kms"^ for WR112, a possible WR+O binary. The high 
velocity for a binary argues for the highest velocity systems 
arising from a combination of DES and BSS. 

A qualitative method to determine the WR runaway 
velocity distribution is to consider the orbital velocities of 
WR bina ries with massiv e O sta r companions. Using the cat- 
alogue of Ivan der Huchtl I 200 ih which lists the masses and 
periods of many WR binaries we find that, without black 
hole kicks, many of these binaries are difficult to unbind be- 
cause the O star is the more massive of the two stars. From 
the binaries we find typically low velocities, similar to those 
predicted in Fig. 3, Oof 10 to 50kms-i for a 3OM0 com- 
panion and up to 70kms~^ for a 10A/q companion. Again 
this leads us to assume that fast and massive WR runaways 
are most likely to be the result of DES. 

3.1.2 Red supergiant runaways 

Runaway red supergiants have a similar distribution to the 
OB runaway stars. It is an interesting question whether such 
stars could be observed and what would be the structural ef- 
fects on their envelopes. The most obvious difference may be 
the distribution of the mass loss in a trail as the star moves in 
space. A trail has recently been discovered for the red giant 
Mira which is a low-to intermediate-mass star with a space 
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Figure 4. The distribution of runaway velocities for neutron 
stars. The smooth black line is the input distribution of neutron 
star kicks. The lines have been normalised so the total probability 
is 1. 
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SN. Their velocities reflect the orbital velocity of the neutron 
star in a binary. Thus, when the neutron star kick distribu- 
tion is estimated from observations, account must be taken 
of the fact that some of the observed neutron star space ve- 
locities are not determined by the neutron star kick. If this 
is not considered the mean kick velocity may be underesti- 
mated and the variance overestimated. Also the upper ve- 
locity side of the population is lower than expected. This is 
because when a binary is unbound as the components sepa- 
rate their mutual gravitational attraction must be overcome. 
Therefore the neutron stars appear to be travelling slower 
than they were at birth. We flnd that we would have to in- 
crease the root mean squared velocity used in the Maxwell- 
Distribution from 265km s~^ to 275km s~^. We note that 
here, we do not consider the possibility of a population of 
low-kick neutron stars emerging from electron-capture su- 
pernovae, as advocated by [Podsiadlowski ct al. (2004). 

Neutron stars are our most common runaway compact 
remnant. Black holes are less populous. We see from Fig. [5] 
that black holes have a distribution with a peak at around 
60 km s~^ and most black holes below 250 km s~^ . Our rarest 
runaway remnants are white dwarfs. The distribution is sim- 
ilar to that for all stars in Fig. [51 This is because white dwarfs 
do not receive any kick and their velocity distribution is rep- 
resentative of their initial binary separations. 

We note that with a typical galactic escape velocity of 
a few times lOOkms"^ we can conclude that most white 
dwarfs and black holes are retained within their host galax- 
ies, while neutron stars are most likely to escape to the in- 
tergalactic medium. 



3.2 Runaway binaries 

The next objects to consider are the runaway binary systems 
that must contain one or more compact objects. Only 20 per 
cent of binaries survive the first supernova. Typically these 
objects have lower velocities than single stars because any 
kick that did not disrupt the binary must be small and had 
to transfer momentum to the system rather than just the 
neutron star. 



Figure 5. Similar to Fig. |4]but for black holes and white dwarfs. 



velocity of 130 kms"^ (|Martin et al.ll2007l 'l. Also the closest 
red supergiant to the Sun, a Ori or Betelgeuse, is a run- 
away star, which mov es with a velocity of about 30kms^^ 
through its local ISM (lUeta et al-lbOOSl ). 

3.1.3 White dwarfs, neutron stars and black holes 

Finally we consider the velocity distributions for single com- 
pact remnants (Fig. [3] and [S]). Their relative formation rates 
and their average space velocities are shown in Tables [S] [71 
[8] and [9] We see that it is most likely that stellar remnants 
are single objects. 

The neutron star velocity distribution is dominated by 
the kick velocity distribution we assume for these objects. At 
low velocities there is a small excess of neutron stars com- 
pared to the input kick velocity distribution. These stars 
come from binary systems that are not disrupted in the first 



3.2.1 Stars with a compact companion 

If a binary survives the first SN then its velocity relative 
to the pre-SN centre of mass changes. The neutron star or 
black hole kick determines this velocity. We see our results in 
Fig. [6] and there are clear differences between the two types 
of remnant. The larger neutron star kicks lead to greater 
velocities so the peak velocity is slightly higher, 60kms~^ 
rather than 20 kms~^ for black hole binaries. The maximum 
velocities possible by each binary reflect this trend with most 
black hole binaries having velocities below 200kms~^ while 
neutron star binaries reaching velocities up to 300 km s^^. 

For both compact remnants a reasonable number of 
these binaries have velocities over 30kms~^. Therefore a 
prediction of our code is that such binaries are less likely to 
be observed within their natal clusters and a r e more likely 
to be runaway stars (e.g. [Mirabel et al. I [2004 [Fragos et al.l 

iooi). 
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Figure 6. Similar to Fig.|4]but for binaries that contain either a 
neutron star or black hole. The lifetime of these systems is taken 
into consideration when calculating the distribution. 
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Figure 7. Similar to Fig. [4] but for systems containing two com- 
pact objects. 



3.2.2 Double compact systems 

Systems which experience two SNe are very unlikely to re- 
main bound. There is very little mass left in the system at 
the time of the second SN, so the system is easily unbound. 
However the more massive the remnant from the first SN the 
more likely the binary is to remain bound but the slower the 
resultant binary will move. The relative birth rates of differ- 
ent single and binary compact remnants are listed in Tables 
[6] and [8] and this trend is clear with there being more black- 
hole primary compact binaries than neutron-star primary 
binaries. Furthermore in Tables [7| and [5] the mean velocities 
of neutron-star primary systems is greater than that when 
the primary produces a black-hole. 

Further insight into the nature of compact binaries can 
be seen in Fig. [T] There is a trend through the distributions 
from only low space velocities for double black-hole binaries 
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Figure 8. Similar to Fig. |4]but for systems containing a compact 
object and a white dwarf. 



to the highest velocities for double neutron-star binaries. 
This demonstrates the most important fact for the remnant 
velocities is the number of large neutron-star kicks a binary 



As mentioned in Section 13.1.31 these velocity distribu- 
tions imply that double neutron-star compact binaries may 
be able to escape their host galaxy because of they experi- 
ence at least one large kick compared to double black hole 
binaries which should be retained by most galaxies. 



3. 2. 3 White dwarf systems 

For completeness we have also included the systems that 
contain a white dwarf and a neutron star or black hole. We 
find in most of these systems the white dwarf forms second. 
The velocity distribution is determined primarily by the kick 
given to the neutron star or black hole at its formation. 

Systems that form white dwarfs first and then neutron 
stars or black holes due to mass transfer are rare. This is to 
be expected because white dwarfs have low mass and, unless 
the kick is small or occurs in a small restricted set of direc- 
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tions, the system becomes unbound. The runaway velocity 
of such systems is greater than systems that experience the 
SN before the WD is formed. This is because while the sys- 
tems are easier to unbind because of their low system mass 
they are also easier to accelerate to high velocities. Therefore 
the peak velocity of these systems is higher and provides an 
observational signature of such systems. The fact that we 
model some such systems indicates that we are producing 
SNe from secondary stars that accrete material to end their 
life in SNe rather than as white dwarfs. 



3.3 Supernovae and gamma-ray bursts 

With the velocities of the runaway stars and their lifetimes 
as runaways we calculate the distance that stars may travel 
before their explosive deaths. This may help to qualitatively 
interpret recent observations concerning the environments 
and locations within the host-galax y of core-collapse SN e 
and gamma-ray bu r sts. F or example iFruchter et al. ll2006ll : 
I James fc Andersonl lliooi): iKellv. Kirshner fc Pahrd (|2008l ') 



and Anderson fc Jamed (|2008l l resolve the host galaxy at 
resolutions of hundreds of parsec. Determining how differ- 
ent superno vae types are d i stribu ted relative to the light 
of a galaxy. iHammer^ et all (|2006l ) studies the host galax- 
ies with greater resolution, calculating the distance between 
long-GRBs and the nearest region containing WR stars. We 
calculate the distance a progenitor travels by multiplying 
the runaway velocity of the progenitor by the time spent as 
a runaway. To calculate the apparent distance observed on 
the sky we also take account of the random orientation of 
the plane of the binary and the phase of the two stars in 
the orbit. When we calculate the apparent distance we also 
take into consideration the three dimensional motion of the 
runaway rel ative to the original binary orbit as given by the 
formulae of lTauris fc"T akcns (1998). 

We analyze our results in two different forms. First we 
arrange the SN progenitors by the type of SN we predict the 
star produces and secondly we consider how distance varies 
with initial mass. Our results are shown in Tables [10] and 
I 111 showing the predicted ratios, mean effective initial mass, 
runaway velocity, distance travelled and ages before core- 
collapse for the different SN types and progenitor sources. In 
calculating the distance travelled we neglect the galaxy po- 
tential. Therefore the furthest dist ances that we predict may 
be modified if this was considered (jVoss fc Taurisll2003l l. We 
note the distance distributions are non-Gaussian and thus 
the mean distances give here should be considered with care. 
In these tables there are also some surprises. For example 
in Table [TT1 the mean mass of type Ic SNe decreases relative 
to the higher metallicity mass in Table 1101 This is primar- 
ily due to the inclusion of QHE changing the outcome of 
lower mass stars and also that the highest mass stars retain 
enough helium to become type lb rather than type Ic (see 
Figure 1). 

Using the method outlined in Section 2.4 above we are 
able to create Fig. [9l and 1 101 showing the distances travelled 
by progenitors of different SN types from our binary mod- 
els alone. There are also panels showing how the different 
progenitor masses contribute to the combined distributions. 
We only consider binaries because we have not included DBS 
runaways in our model. If the number of runaways is similar 
from DES as for BSS then the distribution of SNe for the 



latter should be similar to that for the combined runaway 
population. 

On first inspection of Figures [9] and [10] we see that most 
SNe occur at their initial location. This is no surprise be- 
cause 60 per cent of SNe from binaries are from the primary 
star, with only 40 per cent coming from secondary stars 
which are the most likely to travel any distance from their 
initial location. Primary stars dominate the statistics be- 
cause this also includes merged systems. These increase the 
number of primary stars that have an effective initial mass 
great enough to produce a SN and decrease the number of 
secondary stars available to produce a supernova. The lat- 
ter effect further boosts the apparent contribution of single 
stars to the total number of SNe. 

The property that most affects the range of runaway 
distances for a specific SN type is the mass range of the 
progenitors. Type IIP come from stars with masses below 
2QMq, which travel great distances due to their long life- 
time. The other SN types typically have higher maximum 
initial masses and wider mass ranges (see Fig. [TJ which lead 
to the different distributions to those in the Figs. In the so- 
lar metallicity plot in Fig. [9] we see that type lb SNe have a 
similar distribution to type IIP SNe because most progeni- 
tors come from the same initial mass range. This similarity is 
weaker at the lower metallicity because it becomes swamped 
by the stars with QHE. This is also evident from the mean 
initial mass of the different progenitors in Tables [TO] and [TTI 

Type Ic SNe typically occur closer to their initial po- 
sitions than type IIP and lb SNe. It is difficult to find ob- 
servations with which these predictions can be tested. Use 
of the observations of discussed above would require some 
knowledge of the luminosity of the source stellar population, 
star formation history and the galactic potentials through 
which the runaway stars travel. Also, as yet, these studies 
do not consider the selection effects of detecting supernovae 
in luminous and dusty regions of galaxies. For example, type 
Ib/c supernovae are generally more luminous than type IIP 
supernovae, although there is signifi cant diversity in the SN 
population ([Richardson et aLll2002h . If this difference was 
true then type Ib/c would be more easily discovered in lu- 
minous regions of galaxies where type IIPs might be missed. 
Thus some of the observed difference in the distributions 
of different supernova types could be due to the intrinsic 
luminosities of these objects and selection effects. 

When we conside r that the resolution which 
iKellv. Kirshner fc Pahr j l|2008l ) resolves typical galax- 
ies in their samples is only a few hundred parsecs. Therefore 
it is the relative distributions beyond these distances are 
important. Our results in Figures [9]andll0 [ do therefore qual- 
itative ly agree with the findings of iKellv. Kirshner fc Pahrg 
( j2008l ) that, in general terms, type IIP and lb SNe should 
be less correlated with star-formation in a galaxy than type 
Ic SNe, especially in environments with solar metallicities. 
This is because type II and lb SNe travel similar distances 
because the age and mass ranges at which these SNe 
occur overlap. In comparison the ages for type Ic SNe are 
considerably lower. The type Ic SNe do not travel so far 
from their birthplaces because they arise from the most 
massive progenitors. This is in addition to the lifetimes of 
the progenitors of typical IIP and lb SNe being similar and 
longer than those of type Ic SNe. 
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3.4 Long gamma-ray bursts 

In conjunction with studying SNe we consider the predicted 
distribution for long-duration GRBs and we relate this to 
quas i-chcmically homogeneous evolution (Sect. 2.2). While 
lYoon. Langer fc Normag (,2006i ) relate long-duration GRBs 
to the quasi-chemically homogeneous evolution of rapidly- 
rotating single stars, we consider only the binary-induced 
quasi-chemically hom ogeneous evolution as explored by 
ICantiello et"all (|2007l ). This has the advantage that the pre- 
dictions do not depend on the initial distribution function 
of stellar rotation rates. This also means that we predict 
a wi der range of initial ma sses for long-GRB progenitors 
than lYoon. Langer fc Norma n (2006). For their single star 
models stars would lose significant angular momentum be- 
fore core-collapse, even at lower metallicities than we study 
here. In our model because the stars are spun up they can 
retain more angular momentum until core-collapse. Because 
of this we may be overestimating the number of very mas- 
sive (greater than QQMq) long-GRB progenitors. If we have 
fewer of these objects it would increase the mean distance 
travelled by the progenitor stars given in Table [TT] 

We find that long-duration GRBs have a bimodal dis- 
tribution of distances. First, 20 per cent explode without 
travelling any distance. These are from binaries where the 
primary star does not explode or the supernova order has 
been reversed. Then a second distribution travels on av- 
erage a few hundred parsecs before exploding. This sec- 
ond group of long- GRBs are constrained to greater appar- 
ent distances because of the requirement that the binary 
orbit must be perpendicular to the line of sight for the 
GRB to be observed. Their dist ribution is also similar to 
that for type Ic, as suggested bv lFruchter et al.l l|2006l ) and 
iKellv. Kirshner fc Pahr3 (|2008l '). While at the same time as 
shown in Figure [10] the maximum distances possible are 
smaller than those of other SNe. In addition many of our 
possible GRB progenitors travel a few hundred parsecs from 
their initi al location. Th i s is in agreement with the obser- 
vations of iHammer et all (|2006l ) who find such distances for 
some near by GRB progenitors. Initia lly this result may ap - 
pear at odds with the observations of iFruchter et all (|2006l ) 
but as mentioned the maximum distance that can be trav- 
elled would be unresolved in their observations so GRBs 
would be expected to occur associated with regions of the 
most intense star formation. 

Also we note from Fig. [10] that lower mass long-GRB 
progenitors are more likely to travel further than high 
mass progenit o rs. T he nearby long- GRBs considered by 
[Hammer et al.l ([20061 ) ar e relatively w eak compared to nor- 
mal cosmological GRBs ([Wooslev fc Bl o om. 2006 ). It is pos- 
sible that the nearby bursts come from lower mass progeni- 
tors than the high redshift long-GRBs. However such a de- 
duction is uncertain and cannot be investigated with our 
simple model. 

We have also varied our requirement that the final CO 
core mass must be greater than 7Mq for a long-GRB to 
occur. For example if we restrict long-GRB to have CO 
core masses at core-collapse of SM© rather than IMq as 
in our fiducial model then the mean distance travelled by 
long-GRBs increased to to 240 ± 380pc. The mean effec- 
tive initial mass and lifetimes also change to 21 ± 14Mq 
43 ± 25Myrs. These differences are because the long-GRB 



population is more dominated by the lower mass QHE pro- 
genitors. We find qualitatively that if such a population rep- 
resented the observed population then long-GRBs would be 
the most widely distributed explosive transient of all core- 
collapse events. Therefore out predicted long-GRB progeni- 
tor population is strongly dependent on the approximations 

included in our simulation. 

Two studies which used the results of [Fruchter et al.[ 

(|2006l ) t o determine a l imit o n th e progenitors of long-GRBs 
were bv [Larsson et al.l (|2007[ ) and [Raskin et al. ( 2008 ). They 
found that they could explain the results of iFruchter et al.[ 
(2006) if the initial mass of GRB progenitors is greater than 
2OM0 and 43Mq respectively (with lifetimes of lOMyrs and 
5Myrs). They did not consider the effect of runaway stars in 
this analysis. Our results on GRB progenitors in Table [TT] 
show that GRB progenitors may have effective initial masses 
of 41 ± 17Mq (or an initial mass of 27 ± I6M0 before mass 
accretion). Because the stars can accrete mass from their 
primary companion the lifetime of these possible progeni- 
tors is 15 ± 5Myrs, three times the age normally expected 
for stars of this mass. This i s greater than the maximum age 
predicted by Larsson et al.l ([2007) and Raskin et al. (200^ 
for the progenitors of long-GRBs. The greater age is due to 
the QHE increasing the main-sequence lifetime and rejuve- 
n ation of the sta r s. Ou r result therefore agrees with the mass 
of [Raskin et al.[ (|2008[ ) bu t disagrees with thei r age as it is 
too low. While the mass of [Larsson et al.[ (|2007l ) agrees with 
the initial mass of our progenitors before mass accretion is 
considered and is closer to our age estimate. 

Qualitatively, because of the uncertainties of what is re- 
quired to cause a long-GRB and the effect of DES runaways 
on our results, our model is consistent with the observations 
of the locations of SNe and long-GRBs. Quantitative anal- 
ysis however require involving our runaway models with a 
more detailed galactic model that also considers the stellar 
populations (that consider binaries) that give rise to these 
events. 



3.5 Merger events of compact-object binaries 

We have also calculated similar plots as in the previ- 
ous subsection for mergers of the compact remnant bina- 
ries we discuss in Section 13.2.21 We use the general rel- 
ativistic weak- field a ppro ximation equation employed by 
[Hurlev. Tout fc PoQ (|2002i ) to calculate the time for such 
binaries to merge via gravitational radiation. We then only 
count it as a merger if the time we estimate is less than 13.5 
Gyrs. 

Our calculations are based on the present day popula- 
tion. I.e., they do not include a realistic simulation of the 
changing populations over cosmological times. Our results 
should thus hold as long as most systems merge within the 
age of the Universe. As shown below, this is mostly the case. 

We show the distribution of merger times in Fig. 1111 
We see that it is more probable that a system takes a long 
time to merge. However there are also systems that have 
extremely short merger times and so could mer ge within star 
formati on regions. The s e resu lt s mirror those of Dewi fc Fold 
(i2003l). [Voss fc TaurisI ([20q3[ ). [Belczvnski et al.l (|2006l ) and 
[O'Shaughnessv et al.[ (|2008l ). 

Over a similar range of merger times [Belczvnski et al] 

(l 2006l ) have two sharply defined peaks in the merger rate at 
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Figure 9. The distribution of distances travelled by different SN progenitors arranged by SN type and initial mass from our binary 
models alone. The lines and contours are normalised so the total number of SNe is 1. The plots are for solar metallicity. Contours are 
logj^g(probability/MQ pc) 



around 10^ and 10^ years. They found that the two peaks 
correspond to two different compact binary formation sce- 
narios. In the first each star experiences one mass transfer 
episode while in the second the secondary experiences an ex- 
tra mass transfer episode when it becomes a helium giant. 
This second scenario leads to the tightest binaries and thus 
the shortest merger times. We do not find such a sharply 
defined peaks however our range of possible times is compa- 
rable. We may be incorrectly estimating some of the merger 
times due to our basic treatment of the gravitational radia- 
tion and assumption of circular orbits in the stellar models. 

The merger rates, shown in Table [S] and Fig. 1111 in- 
dicates that merging black holes and neutron stars are the 
most common visible mergers. However we are not able to es- 
timate the merger rates for systems containing white dwarfs 
due to limitations of our grid of models for the secondaries 
evolution. We are able to calculate the birthrate and give 
this in Table [S] If a similar number of these systems merge 
within the age of the universe as for the other compact bi- 
naries then the rate of white dwarf-neutron star or black 



Table 5. Predicted rates of binary compact object mergers and 
births relative to the rate of core collapse supernovae, at the two 
considered metallicities. 



Z WD-NS 


WD-BH 


NS-NS 


NS-BH 


BH-BH 


Compact-system 


merger rate 








0.004 




0.00007 


0.0017 


0.0019 


0.020 




0.00009 


0.0008 


0.0021 


Compact-system 
0.004 0.0045 
0.020 0.0054 


birth rate 
0.0078 
0.0074 


0.0003 
0.0004 


0.0052 
0.0029 


0.0292 
0.0181 



hole mergers may be greatest. The rate of such mergers 
must be less than l/80t h of the core collapse supernova rate. 
iThompson et all (|20o3) suggest a rate between l/20th and 
l/40th of the SN rate from the observation of one such sys- 
tem. Because the observable outcome of such events is not 
yet well studied l|Paschalidis et al]|2009l ). it seems unclear 
whether they have observable counterparts. We suggest that 
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Figure 10. Similar to Figure [9] but now SMC-like metallicity. Tlie red contours in tlie lower right panel indicate that distribution for 
long-GRB progenitors. 
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Figure 11. The distribution of time required for compact binaries to merge. The lines are normalised so the total number of mergers is 
one. 
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it is possible for these systems the merger is most likely to be 
an extended period of mass transfer. When the white dwarf 
is less massive th an the neutron st ar any mass transfer is 
sufficiently stable (|Tout et al.ll 19971 ) so these binaries might 
be observed as X-ray binaries with mass-transfer driven by 
gravitational radiation. 

When we consider where compact binary merger events 
might be observed we should bear in mind the results dis- 
played in Fig. [71 Here we find that many compact bina- 
ries have high space velocities and can escape their host 
galaxies, especially double neutron-star systems. Therefore 
we should expect some compact-object binary merger events 
to occur some distance away from any nearby galaxy. There 
is growing evidence that som e such events do occur (jLevanl 
l2007l : ISalvaterra et aLlboioD . Merging double-neutron star 
systems as well as neutron star-black hole mergers are pos- 
sibly related to short gamma-ray bursts. In this case these 
rates should be considered upper limits on the observable 
rate because the progenitor rotation axis must be within a 
small angle to the line of sight to be observed, just as for 
the long-GRBs. 



4 DISCUSSION & CONCLUSIONS 

Our population synthesis of runaway stars is able to predict 
a velocity distribution similar to that observed for OB stars. 
We predict that in the Galaxy the number of all O stars that 
are runaways fro m ESS is between . 5 to 2 .2 per cent. Using 
the catalogue of iMaz-Apellniz et al.l (|2004| ) we estimate the 
observed number is 4 ± 1 per cent, however both numbers 
are uncertain and very sensitive to the total number of O 
stars in the sample and how runaways are identified. 

While we do not consider dynamical ejection from clus- 
ters here, we suggest that this distribution may be similar 
to that predicted for the binary supernova ejection scenario. 
This is because stars ejected by dynamical interactions in 
a cluster are normally the result of binary interactions so 
that the determining factor for both scenarios is the initial 
separation distribution of binary stars. The initial binary 
period/separation distribution has a direct effect on the ve- 
locities of runaway stars. A wide range of binary separations 
is required to explain the observed velocity distribution. A 
distribution that is flat in the logarithm of the separation re- 
produces the observed range of runaway velocities. However 
this separation distribution may underestimate the number 
of close binaries. 

To match the fact that there are a number of observed 
WR runaways with velocities greater than 30kms~^ we in- 
clude black hole kicks when black holes are formed by core- 
collapse. This has implications for the velocities of binaries 
that include a black hole. These should be more often ob- 
served as runaways. For example the X-ray binary Nova Sco 
which has a space velocity of 150 ± 19kms~'^ provides evi- 
dence that at least some black hol es may have kicks at their 
formation. iNelemans erahl (|l999l ) suggest that such large 
velocities for black-hole binaries do not require a large kick. 
From the black-hole binaries they list we find the mean ve- 
locity is 34 ± 32km s^^, this is similar to the distribution of 
synthetic black hole binaries shown in Figure [S] 

We consider the final outcome of our stars and the ve- 
locity distribution of compact stellar remnants. Our models 



predict that there should be a small surplus of neutron stars 
with space velocities less than 200kms~^. These require a 
system not to be unbound by the first SN but by the second. 
If these were not taken into account when the neutron star 
kick is determined from observations it would lead to an un- 
derestimate of the kick velocity. Most white dwarfs and black 
hole runaways have velocities typically below those expected 
for neutron stars. The final outcome of evolution of all our bi- 
naries is that most compact stellar remnants are single. Only 
20 per cent of all binaries that experience at least one SN 
remain bound. The rarest are white dwarf and neutron star 
binaries where the white dwarf formed first. Normally such 
systems can be identified by their highly eccentr ic orbits, for 
exam nle PSR B2303-H46 and PS R J1141-6545 (|Kaspi et all 
I2OOOI : Ivan Kerkwiik fc Kulkarni|[l999h . There is a predicted 
difference in the space velocity between such binaries and 
those in which the neutron star forms first. 

Whether a SN progenitor is a runaway star is an im- 
portant factor in where a SN might occur in a galaxy. We 
have shown that less massive stars travel further from their 
initial positions than more massive stars. However the ma- 
jority of stars explode without travelling a large distance. 
When we consider the different SN types we find that be- 
cause type II and type lb SNe come from progenitors less 
massive than 20Mq they can travel furthest from their ini- 
tial positions, 48 ± 210pc and 34 ± 150pc respectively on 
average. While type Ic progenitors come from more massive 
stars and therefore have shorter lifetimes and travel shorter 
distances, 17±74pc. This picture changes at low metallicity 
with the inclusion of quasi-homogeneous evolution, which 
prolongs the life and the distance to which type Ic and type 
lb progenitors can travel. 

Nearly one third of long-GRB progenitors may travel 
a few hundred pc from their birthplaces before experienc- 
ing core-collapse. But, they can travel shorter distances 
and/or explode before type IIP and lb SNe and so should be 
the most tightly correlated with star-formation in the host 
galaxies. These distances however are dependent on the pa- 
rameters taken for a model to give rise to a long-GRB. 

Type IIP supernova progenitors that are secondaries in 
a binary are distributed throughout a galaxy due to BSS 
more than progenitors of long-GRBs. These remain closer 
to their birth places and have shorter lifetimes so explode 
when the star- f ormat i on episode the y belong to is still visible 
iLarsson et"al] (|2007^ : [Raskin et all ((200^). Therefore these 
effects in combination can explain the observed distribution 
of SN within galaxies. To infer the nature of long-GRB pro- 
genitors from such observations runaway stars must be con- 
sidered. We find that of the supernovae in Tables [TOI and [TTl 
20 per cent of the SN come from BSS runaway stars. If there 
are at least as many DES runaway progenitors then between 
20 to 40 per cent of SNe has a progenitor that was runaway 
star. A similar fraction w as suggested from observations by 
I James fc AndersonI (| 2006V Without considering the effect of 
runaway stars, the mass of long-GRB progenitors is likely to 
be overestimated. 

Neutron star-black hole binaries and more so double 
neutron star systems, which are candidate short-GRB pro- 
genitors, are found to reach large runaway velocities. We 
estimate their merger timescales and find that they should 
often merge outside of their host galaxy. Compact mergers 
of a white dwarf with either a neutron star or black hole 
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are estimated to be at least as common as other compaet 
binary mergers. Due to their low runaway velocities and 
merger timescales, they should mostly merge within their 
host galaxy. The observational signature of such events ap- 
pears to be as yet unclear. 
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Table 6. The relative birth rates of binary- produced unbound compact objects (white dwarfs, neutron stars and black holes), and 
compact objects which remain part of a binary system, compared to the birth rate of single compact objects (last line). For compact 
objects remaining in binaries, it is distinguished whether they form from the primary or the secondary star. These results refer to our 
population synthesis model at a metallicity of Z=0.020. The columns indicate the remnant from the primary star of the binary with the 
rows referring to the remnant from the secondary star. 





Unbound 


Primary WD 


Primary NS 


Primary BH 


Unbound 




0.0176 


0.1397 


0.0747 


Secondary WD 


0.0280 


0.1634 


0.0032 


0.0045 


Secondary NS 


0.0916 


0.0002 


0.0003 


0.0013 


Secondary BH 


0.0354 


0.0001 


0.0005 


0.0112 


Single 




0.1831 


0.1481 


0.0972 



Table 7. The mean space velocities of different single and binary compact remnants at a metallicity of Z=0.020. 





Unbound 


Primary WD 


Primary NS 


Primary BH 


Unbound 




11 ± 3 


420 ± 180 


100 ± 70 


Secondary WD 


24 ± 18 





81 ±38 


51 ± 27 


Secondary NS 


390 ± 180 


140 ± 55 


160 ± 62 


110 ±45 


Secondary BH 


100 ± 67 


120 ±53 


110 ±52 


53 ±32 


Single 







430 ± 190 


96 ±61 



Table 8. The relative birth rates of different single and binary compact remnants at a metallicity of Z=0.004. 





Unbound 


Primary WD 


Primary NS 


Primary BH 


Unbound 




0.0230 


0.1681 


0.0763 


Secondary WD 


0.0220 


0.1175 


0.0032 


0.0048 


Secondary NS 


0.0897 


0.00005 


0.0002 


0.0024 


Secondary BH 


0.0461 


0.0010 


0.0014 


0.0214 


Single 




0.1279 


0.1887 


0.1063 



Table 9. The mean space velocities of different single and binary compact remnants at a metallicity of Z=0.004. 



Unbound Primary WD Primary NS Primary BH 



Unbound 








10 ±2 


420 ± 180 


95 ± 


72 


Secondary WD 


26 


± 


20 





81 ±39 


56 ± 


30 


Secondary NS 


380 


± 


180 


89 ± 29 


160 ± 65 


97 ± 


46 


Secondary BH 


93 


± 


55 


76 ± 36 


87 ±55 


37 ± 


28 


Single 











430 ± 180 


94 ± 


71 
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Table 10. The relative rates of the different SN types and their progenitors mean effective initial mass, mean space velocity and mean 
distance travelled from point of origin at a metallicity of Z=0.020. The relative rates and means are also split up into the values for 
single stars, primaries and so( ondarics. 







IIP 


non-IIP 


lb 


Ic 


L-GRB 


Fraction 


All 


0.582 


0.119 


0.068 


0.231 


0.0 




Single 


0.287 


0.046 


0.007 


0.058 


0.0 




Primary 


0.150 


0.055 


0.045 


0.123 


0.0 




oc^^ii^cLi y 


0.144 


0.018 


0.016 


0.050 


0.0 


< Mi > 


All 


10.7 ±3.0 


17.7 ±6.3 


16.8 ±5.8 


37.1 ± 22.8 






Single 


11.1 ± 2.9 


21.9 ±2.3 


27.5 ±0.5 


50.7 ±21.8 






Primary 


10.1 ±2.9 


13.2 ±5.4 


15.6 ±4.5 


32.7 ±22.5 






Secondary 


10.5 ± 3.2 


20.3 ±6.6 


15.6 ±5.8 


32.3 ± 18.0 




< f run > 


All 


4.4 ± 12.0 


2.0 ±6.1 


5.0 ± 13.4 


4.0 ± 10.8 




[kms-1] 


Single 


0.0 


0.0 


0.0 


0.0 






Primary 


0.2 ± 1.2 


1.0 ±3.5 


0.3 ± 1.8 


0.4 ±1.9 






Secondary 


17.5 ± 18.7 


10.0 ±11.6 


20.3 ±21.0 


17.4 ± 17.4 




< drun > 


All 


48 ± 210 


10 ±39 


35 ± 150 


17 ±74 




[pc] 


Single 


0.0 


0.0 


0.0 


0.0 






Primary 


2 ± 16 


8 ±36 


2± 15 


1 ± 10 






Secondary 


190 ± 380 


38 ± 71 


150 ± 280 


74 ± 150 




^ trun ^ 


All 


36 ±31 


22 ± 38 


21 ±21 


9 ± 7 




[Myr] 


Single 


26± 10 


9±1 


7.0 ±0.1 


5±1 






Primary 


33 ± 15 


26 ± 20 


17 ±9 


10 ± 5 








.^)7 ± r,2 


1() ± 85 


:i7 = :-! 1 


l:-! ± 11 





Table 11. The relative rates of the different SN types and their progenitors mean effective initial mass, mean space velocity and mean 
distance travelled from jjoint of origin at a metallicity of Z=0.004. The relative rates and means are also split up into the values for single 
stars, primaries and secondaries. The mean mass in brackets for the long-GRBs is the actual initial mass before the stars have accreted 
material from the primary. 







IIP 


non-IIP 


lb 


Ic 


L-GRB 


Fraction 


All 


0.605 


0.141 


0.164 


0.091 


0.001 




Single 


0.324 


0.055 


0.025 


0.000 


0.0 




Primary 


0.163 


0.079 


0.079 


0.055 


0.0 




Secondary 


0.118 


0.008 


0.060 


0.036 


0.001 


< Mi > 


All 


9.9 ±3.4 


21.4 ± 10.0 


30.0 ±25.3 


25.6 ± 18.5 


40.7 ± 16.8(15.2 ± 12.1) 


[Mq] 


Single 


10.6 ±3.5 


28.4 ± 5.9 


66.2 ±20.4 








Primary 


9.2 ±3.1 


16.3 ±9.5 


28.5 ±24.0 


22.8 ± 17.9 






Secondary 


9.0 ±2.8 


24.4 ±5.4 


17.3 ± 11.1 


30.0 ± 18.4 


40.7 ±16.8(15.2 ±12.1) 


< Vrun > 


All 


4.2 ±12.1 


1.0 ±3.6 


6.0 ± 14.0 


8.6 ± 16.3 


12.0 ± 10.9 


[kms-1] 


Single 


0.0 


0.0 


0.0 








Primary 


0.1 ±0.9 


0.6 ±2.1 


0.3 ±1.6 


0.2 ± 1.2 






Secondary 


21.2 ± 19.9 


12.7 ±6.5 


15.8 ± 19.4 


21.3 ±19.9 


12.0 ± 10.9 


< rfrun > 


All 


54 ± 240 


3 ± 19 


75 ± 250 


60 ± 150 


88 ± 100 


[pc] 


Single 


0.0 


0.0 


0.0 


0.0 






Primary 


0.5 ±4 


3 ±21 


3 ±21 


1±6 






Secondary 


280 ± 480 


27 ±35 


200 ± 390 


150 ± 210 


88 ± 100 


*C trun ^ 


All 


39 ±25 


15 ±12 


27 ±30 


16 ±7 


15 ±5 


[Myr] 


Single 


31 ±14 


8±1 


4±1 








Primary 


41 ±20 


21 ±14 


17 ±15 


15 ±7 






Secondary 


58 ±40 


9±2 


49 ±36 


17 ±7 


15 ±5 



20 John J. Eldridge 



Vink J.S., de Koter A., Lamers H.J.G.L.M., 2001, A&A, 
369, 574 

Voss R., Tauris T.M., 2003, MNRAS, 342, 1169V 

Wellstein S., Langer N., 1999, A&A, 350, 148 

Woosley S.E., Bloom J.S., 2006, ARA&A, 44, 507 

Yoon S.-C, Langer N., 2005, A&A, 443, 643Y 

Yoon S.-C, Langer N., Norman C, 2006, A&A, 460, 199Y 

Yoon S.-C, Woosley S.E., Langer N., 2010, 

arXivl004.0843Y 
Zinnecker H., Yorke H.W., 2007, ARA&A, 45, 481Z 



6 APPENDIX 

In this appendix, we present additional results from our sim- 
ulations. These are placed here not to divert the main results 
we wished to cover. These Figures contain similar informa- 
tion to that presented in the main paper. However these pro- 
vide greater clarity on some specific results of the simulation. 
Fig. [12] is similar to Fig. [3] but with the range of velocities 
included widened to cover velocities less than those typically 
expected for runaway stars. Fig. [TH] is similar to Fig. O but 
is over a linear distance scale. Fig. [14] is also similar to Fig. 
[9] but is for cummulative probability. Then Figs. [15] and [16] 
are also similar to these plots but split the SN events up into 
different initial-mass ranges rather than SN types. Finally 
Figs. [17] are similar to the other plots but show results for 
the compact-object mergers rather than SNe. 

This paper has been typeset from a T^jX/ I^Tp^X file prepared 
by the author. 
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Figure 12. Similar to Figure [3] but witli the velocity range below the normal velocity limit required for a star to be identified as a 
runaway. 
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Figure 13. Similar to the upper left panels in Figures l9l and [TOl but with a linear scale on the x-axis. 
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Figure 14. Similar to the upper left panels in Figures l9land [T0l but with the cummulative probability with increasing progenitor distance 
on the y-axis with a linear scale. 




Figure 15. Similar to the upper left panels in Figures l9l and [Tol but with a linear scale on the x-axis and the lines now represent different 
initial mass ranges for the progenitor stars. 
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Figure 16. Similar to ttie upper left panels in Figures [9l and [TOl with the cummulative probability with increasing progenitor distance 
on the y-axis with a linear scale with a linear scale. 
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Figure 17. Similar to the upper left panels in Figures [9] but now for the mergers of compact objects shown in Figure [TT] with the 
cummulative probability of a merger versus increasing progenitor distance on the y-axis with a linear scale with a linear scale. 



